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Using helium atom scattering Hulpke and Liidecke recently 
observed a giant phonon anomaly for the hydrogen covered 
W (110) and Mo (110) surfaces. An explanation which is able 
to account for this and other experiments is still lacking. Be- 
low we present density-functional theory calculations of the 
atomic and electronic structure of the clean and hydrogen- 
covered Mo (110) surfaces. For the full adsorbate monolayer 
the calculations provide evidence for a strong Fermi surface 
nesting instability. This explains the observed anomalies and 
resolves the apparent inconsistencies of different experiments. 
PACS numbers: 68.35.-p,73.20.-r,73.20.Mf 

Recently, giant anomalies in the surface phonon spec- 
tra have been detected by hehum atom scattering (HAS) 
on the (110) surfaces of tungsten and molybdenum, when 
these are covered with a full monolayer of hydrogen ^^ ■ 
One critical wave vector Q"^"^ is along the [001] direction 
(PH) and has a length of 0.95 A'^ for W and 0.90 A"! 
for Mo. Anomalies have also been found for directions 
deviating from PH. Along PS they are located at S with 
a critical wave vector Q^^ (IQc2''l — 1.225 A~^ for Mo 
and 1.218 A~^ for W) |j]. The experimental dispersion 
curves (see Ref. H) are only available for PH and PS. 
Thus, we focus our attention particularity on these two 
directions. For both surfaces two different softenings in 
the surface phonon branches were observed at the critical 
wave vectors: a smaller dip and a very deep and sharp 
indentation from Tiuj « 15 meV to fiLO ~ 2 meV. Only the 
phonon dispersion curves of quasi-one-dimensional con- 
ductors, like KCP [K2Pt(CN)4Br3] ||] are characterized 
by similarly deep anomalies strongly localized in recip- 
rocal space. A pronounced but less sharp damping of 
longitudinal surface phonons was observed for the (100) 
surfaces of W |^,0 and Mo [g[ and was explained in terms 
of the nested structure of the Fermi surface [§,0 . In the 
case of W (100) this Fermi surface nesting is apparently 
so strong that it induces a rebonding at the surface and 
a structural rearrangement, i.e. a c(2 x 2) surface re- 
construction at temperatures below Tc « 250 K (for a 
review see |1^). For the clean (110) surfaces of W and 
Mo no indication of an anomaly is found, but this only 
appears when the surfaces are covered with a full layer 
of hydrogen. Nevertheless, a close link between the sur- 
face phonon anomalies and H vibrations seems to be ruled 
out, since the HAS spectra remain practically unchanged 
when deuterium is adsorbed instead of hydrogen P,p| . 

Angular resolved photoemission (ARP) studies have 
been performed by Kevan and coworkers |12-d4| for clean 
and hydrogen covered W(llO) and Mo (110). These 



studies give no evidence of the existence of parallel seg- 
ments of the Fermi-surface contours separated by wave 
vectors comparable with the HAS determined critical 
wave vectors and thus there appeared to be compelling 
reason for abandoning the nesting mechanism as a pos- 
sible origin of the anomalies seen in HAS. 

The situation became even more puzzling when very 
careful high resolution electron energy loss spectroscopy 
(HREELS) studies ||l|l of hydrogen covered W (110) ob- 
served the Rayleigh wave phonon branch with the small 
dip at Qci'' exactly as the HAS results. However, the 
giant anomaly was not detected. 

At present there exists no explanation which is able to 
account for the observed two phonon anomalies and is 
compatible with the different experimental data (HAS, 
HREELS, and ARP). We therefore performed ab-initio 
calculations of the surface atomic geometries and the 
electronic properties of the clean and H-covered Mo (110) 
surface. The results for the adsorbate covered surface 
provide clear evidence for a Fermi surface nesting in- 
stability at wave vectors which are in remarkable agree- 
ment with those at which the phonon anomalies occur 
for this system [Q*'i = (0.86 ± 0.02, 0) k-\ Qlf = 
(0.90,0) A-i, Q;:^ = (1.00 ± 0.02,0.71 ± 0.02) A-\ 
Qlf = (1.00,0.707) A-i]. We choose the a;-axis along 
[001] and the y-axis is along [TlO]. The corresponding di- 
rections in reciprocal space are PH and FN, respectively. 
The calculations were performed using density-functional 
theory together with the local-density approximation for 
the exchange-correlation energy functional jig] • The full- 
potential linearized augmented plane wave method |17| ] 
is employed which we enhanced by the calculation of 
forces |1^. This enables an efficient evaluation of the 
multi-layer substrate relaxation and the atomic positions 
at the surface. The substrate is modeled by a seven layer 
slab repeated periodically with a separation of 8.8 A of 
vacuum. The k-integrations are evaluated on a mesh of 
64 equally spaced points in the surface Brillouin zone 
(SBZ). The mufhn-tin radii are chosen to be 1.27 A and 
0.48 A for Mo and H, respectively. The kinetic-energy 
cutoff for the plane wave basis needed for the interstitial 
region is set to 12 Ry, and the {I, m) representation (in- 
side the muffin tins) is taken up to /max = 8. For the 
potential expansion we use a plane- wave cutoff energy of 
64 Ry and a (l,rn) representation with Zmax = 4. The 
core states as well as the valence states are treated non- 
relativistically. The calculated in-plane lattice constant 
is 3.13 A without including zero point vibrations, which 
compares well with the measured bulk lattice parameter 



(3.148 A at room temperature [Q). The values of the 
surface relaxation for the clean surface are in excellent 
agreement with those calculated by Methfessel et al. pc| ] 
who employed the full-potential linearized muflin-tin or- 
bitals approach. 

For a meaningful theoretical analysis of the surface 
electronic structure it is important to use the correct 
atomic geometry. For the hydrogen adsorption system 
at a full monolayer, i.e. at the coverage for which the 
phonon anomaly was observed experimentally, neither 
the hydrogen adsorption site nor the substrate relaxation 
have been determined so far. In Table I we summarize 
the calculated adsorption energy differences A£^ad as well 
as the relaxation parameters for all possibly important 
hydrogen sites. The geometries are defined in Fig. ||. 
The results identify the hollow site as the most stable 
position. The [110]-offset from the long-bridge position 
is yn = 0.55 ± 0.01 A. Within the computational error 
the hydrogen relaxes exactly in the geometric threefold 
site and there is no theoretical evidence for a pronounced 
top-layer-shift reconstruction in agreement with the low- 
energy electron diffraction experiments pl[ . In compar- 
ison both bridge positions are clearly energetically un- 
favorable, and the on-top position is even worse. The 
theory predicts that the clean-surface inward relaxation 
of —4.5 %do, where do is the bulk inter-layer spacing, 
is reduced for the long-bridge, short-bridge, and hollow 
geometries of the adsorbed hydrogen. In the case of the 
on-top adsorption site we obtain an outward relaxation. 
Note that the numerical error for all calculated inter-layer 
spacings is about ±0.3 %do- 

In Fig. we show the Fermi surfaces for the clean 
Mo (110) surface and for the H-covered surface. On the 
left the theoretical results of the iV-particle ground state 
are plotted and the right side displays the experimental 
ARE data of Jeong et al. |1^ . For the adsorbate system 
we present the results for the hollow position. 

The shaded areas are the projection of the bulk Fermi 
surface onto the (110) SBZ. For the experimental anal- 
ysis it was calculated using a tight-binding interpolation 
scheme [g2[. Differences between this projection and the 
one of the present work may arise from the lack of self- 
consistency of the tight-binding approach. For the clean 
surface (see Fig. ||a) the calculations identify four bands, 
which are highly localized at the surface. There is a band- 
circuit centered at F, one centered at S, and one centered 
at N. Furthermore, there is a band more or less parallel 
to the FN direction. The latter one is of (dsz^-r^, dxz) 
character, and also the other surface features are due 
to Mo d bands. All of these bands have been observed 
experimentally, and for all of them the calculated and 
measured k-dependence is in very good agreement (com- 
pare Figs. 0a and b). Clearly worse agreement is found 
when we compare the results for the adsorbate systems 
(see Fig. |^, lower panels). 

At first we discuss the changes of the theoretical Fermi 



surface induced by the hydrogen adsorption (Fig. ga and 
c), and then we return to the experiments. For the 
present discussion the most important effect is that the 
{d3z2_,.2,dxz) band is shifted away from the FN line. A 
significant fraction is now in the stomach shaped band 
gap. This shift is not a consequence of hybridization be- 
tween the (d3^2_^2, dxz) states and hydrogen orbitals but 
it results because the adsorption of hydrogen changes the 
surface potential and shifts the entire ((i3j,2_^2, d^z) band 
to lower energies. As this band disperses upward when k 
increases away from F, the Fermi energy cuts the shifted 
band at a longer k-vector. The hydrogen induced bond- 
ing states are found at lower energies (about 5 eV below 
the Fermi level), and the antibonding states are at higher 
energies (about 4 eV above the Fermi energy). 

Figure gc reveals that the Fermi-surface contours as- 
sociated to the (d322_^2, cia;2;)-like band run parallel to 
the FN direction and perpendicular to FS in significant 
parts of the SBZ giving rise to a quasi-one-dimensional 
nesting. The k-vectors connecting these bands in differ- 
ent sections of the Brillouin zone are Q*i and Q*2. We 
recall that the highest occupied Kohn-Sham level of the 
self-consistent A^-electron calculation equals the electron 
chemical potential, i.e. the Fermi level. However, the 
Kohn-Sham Fermi surface is (in principle) not an ob- 
servable. The agreement between Fig. ga and b and the 
fact that the wave function character does not change 
significantly upon hydrogen adsorption suggests that the 
difference between the Kohn-Sham Fermi surface and the 
true Born-Oppenheimer Fermi surface is not very impor- 
tant for this system. We therefore conclude that the nest- 
ing instability predicted by the Fermi surface of Fig. 0c 
can be trusted, and due to its one-dimensional nature 
we expect important consequences pJ. Such an insta- 
bility may induce a static distortion provided that the 
electronic energy gain is stronger than the elastic energy 
cost needed to displace the top layer atoms. Because the 
(110) surface is the closest packed and most stable bcc 
surface such a Peierls distortion is somewhat unlikely. 
When a static distortion is not realized, there will be 
strong dynamical consequences (similar to a dynamical 
Jahn- Teller effect in molecules) due to thermally excited 
electron- hole pairs with a wavevector equal to the critical 
wavevectors and their coupling to surface phonons. 

Since for metal surfaces helium atoms scatter at a dis- 
tance of about 3-4 A in front of the surface HAS will 
detect electron charge density oscillations close to the 
Fermi level ^. For the system Mo (110) at the crit- 
ical wave vectors these oscillations are associated with 
electron-hole pair excitations involving the {d^z-2_^2,dxz) 
band. We therefore propose that the two bands seen by 
HAS are due to a hybridization between lattice vibra- 
tions and electron-hole pairs. Thus these studies rep- 
resent an experimental manifestation of a nearly one- 
dimensional Kohn anomaly. While the small dip is due 
to a more phonon like excitation, the sharp and giant 



anomaly is due to a more electron-hole like excitation 
of nested electron-hole pairs. The anomalies persist in 
the phonon dispersion curves along directions deviating 
from FH Pfl. This result is directly understandable in 
terms of the theoretical results of Fig. ||c, because nesting 
is present also for the directions in the surface Brillouin 
zone out of the [001] azimuth. Because of the form of the 
(d3z2_r2, dxz) band Fermi-surface contour, the anomalies 
are particularly strong along FS in agreement with the 
available experimental data H. We expect a weakening 
of the indentations for directions between FH and FS be- 
cause of the less effective nesting. A careful reanalysis of 
the experimental data is in progress [p4| . 

Our interpretation of the HAS results also explains 
why in HREELS only the Rayleigh mode with the small 
dip was observed. This technique couples to vibrations 
of the nuclei and is practically insensitive to electronic 
charge density oscillations |g5[. Because the sharp and 
giant anomaly has only little vibrational character, it is 
practically impossible to excite it by high-energy elec- 
trons. 

The important question remains why the ARP mea- 
surements for the hydrogen covered Mo (110) and the 
theoretical Kohn-Sham Fermi surface are so much dif- 
ferent. The results of Fig. He are due to photoelectrons 
which carry the information about the difference of the 
N particle ground state and the N — 1 particle system 
with a hole at |eF, k). Typically an analysis is done under 
the assumption that the electronic system is in its ground 
state and that vibrational excitations as well as interac- 
tions between electronic and vibrational excitations can 
be ignored. The theoretical results shown in Fig. pic, 
which are obtained under exactly this hypothesis, dis- 
close that for this system the assumption is unfounded. 
In fact, the calculations predict highly localized paral- 
lel bands at the Fermi surface with a one-dimensional 
nesting vector. As a consequence, at Qd and Qc2 many 
electron-hole pairs will be thermally excited, and these 
will couple to phonons. This implies a breakdown of 
the Born-Oppenheimer approximation. While the calcu- 
lations identify a strong nesting, the ARP experiments 
investigate a system which has already reacted to this 
instability. When a static distortion is hindered, ARP 
will measure a state with a significant number of low en- 
ergy excitations, and in particular at the critical wave 
vectors the electron and lattice dynamics cannot be de- 
coupled. In fact, we note that the experimental Fermi 
surface of the H covered surface seems to exhibit a dif- 
ferent translational symmetry than that of the (1x1) 
SBZ (see the band circuits at the S point in Fig. 0d). We 
cannot rule out, however, that this is due to inaccuracies 
in the experimental data or analysis. We are not aware 
that for any other surface such a serious breakdown of 
the Born-Oppenheimer approximation (compare Figs. 0c 
and d) has been observed so far. Therefore, some doubts 
remain at this point if an alternative explanation may 



exist which can also resolve the differences of the experi- 
mental ARP data and the calculated results. Obviously, 
the static distortion of the W(IOO) surface is also in- 
duced by a Kohn anomaly, but after the reconstruction 
has taken place the Born-Oppenheimer approximation is 
again valid. 

In conclusion, the calculated Kohn-Sham Fermi surface 
can explain the HAS and HREELS experiments. It is in- 
teresting that the same band which is responsible for the 
pronounced Kohn anomaly identified above for Mo (110) 
is in fact also responsible for notable features of all group 
VIA materials. In particular we mention the spin-density 
wave in bulk Cr and the phonon anomalies along the [001] 
direction in bulk W and Mo. Also the reconstruction of 
Mo (100) and W(IOO) is induced by the nesting of this 
band. What makes the anomaly at the (110) so partic- 
ularly strong is that the adsorption of hydrogen shifts it 
into a region of k-space where it becomes localized at the 
surface. Thus the band becomes two- and the nesting 
one-dimensional. Neither the hydrogen wave functions 
nor the hydrogen vibration are directly involved in the 
resulting anomalies. The effect of hydrogen is mainly to 
change the potential at the surface and to shift the Mo 
surface states. Thus, it is well possible, that other adsor- 
bates can have a similar effect on the surface vibrations 
and electron-phonon coupling. We also expect that for 
the clean surface a weak signature of the {d2^2_j.2,dxz) 
band should exist at k « 0.6 A"-'^. However, for the clean 
surface the band is not a surface state but a broad and 
less localized surface resonance. The thermal excitations 
of electron-hole pairs at Mo(110):H and the suggested 
breakdown of the Born-Oppenheimer approximation call 
for additional photoemission experiments, in which the 
Fermi surface is carefully studied also at lower tempera- 
tures. In fact, also a detailed analysis of the temperature 
dependence of the small dip seen in HAS and HREELS 
would be most interesting. Furthermore, calculations of 
the electron-phonon coupling, electron-hole excitations 
and frozen phonons at these surfaces would be very help- 
ful, but for transition metal surfaces this is not yet fea- 
sible with the required accuracy. 

We thank J. Liidecke, E. Hulpke, and E. Tosatti for 
stimulating discussions. 
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FIG. 2. Fermi surfaces of clean Mo (110) [upper panels] 
and Mo (110) with a full monolayer of hydrogen [lower pan- 
els] . The experimental results are from [jl3| . Shaded areas are 
the projection of the bulk Fermi surfaces onto the (110) sur- 
face Brillouin zone. The solid (dotted) lines of the theoretical 
results denote surface resonances or surface states which are 
localized by more than 60 % (30 %) in the two top Mo layers. 
Qei and Qc2 are the theoretical critical wavevectors. 
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TABLE L Calculated geometries and adsorption-energy 
differences for clean Mo (110) and for the H-covered surface. 
For the latter the results for the long-bridge, short-bridge, 
hollow, and on-top sites (see Fig. 1) are compiled. dn-Mo 
represents the height of the hydrogen atom above the surface 
while Adij is the percentage change of the inter- layer distance 
between the i-th and the j-th layer with respect to the bulk 
inter-layer spacing do. 



FIG. 1. Surface geometry of Mo (110) with the long-bridge 
(L), short-bridge (S), hollow (H) and on-top (T) sites 
marked, j/h is the [IlO]-offset of the hollow position from 
the long-bridge site. 
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